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Abstract 

Quasielastic neutrino-nucleon scattering is a basic signal process for neutrino oscillation 
studies. At accelerator energies, the corresponding cross section is subject to significant 
uncertainty due to the poorly constrained axial-vector form factor of the nucleon. A 
model-independent description of the axial- vector form factor is presented. Data from the 
MiniBooNE experiment for quasielastic neutrino scattering on 12 C are analyzed under 
the assumption of a definite nuclear model. The value of the axial mass parameter, 
rriA = 0.85j[Q'Q7 ± 0.09 GeV, is found to differ significantly from extractions based on 
traditional form factor models. Implications for future neutrino scattering and pion 
electroproduction measurements are discussed. 



1 Introduction 



High statistics neutrino experiments are probing the hadronic structure of nuclear targets at 
accelerator energies with ever greater precision. Extracting the underlying weak-interaction 
parameters, or new physics signals, requires similar precision in the theoretical description of 
the strong interactions. 

A basic cross section describes the charged- current quasielastic scattering process on the 
neutron, 



Recent evidence indicates a tension between measurements of this process in neutrino scat- 
tering at low [TJ [2J El H] and high [5] neutrino energies, and between results from neutrino 
scattering and results inferred from pion electroproduction [BJ. In particular, with a com- 
monly used dipole ansatz for the axial-vector form factor of the nucleon, 



different experiments have reported values for the so-called axial mass parameter m A ipoe . 
World averages reported by Bernard et al. [6J find comparable values obtained from neutrino 
scattering results prior to 1990, m d f olc = 1.026 ±0.021 GeV, and from pion electroproduction, 
m di P oic = ^ L069 _ Q55 ^ ± o.oi6 GeVQ The NOMAD collaboration reports m^ ipole = 
1.05 ± 0.02 ± 0.06 GeV. In contrast, MiniBooNE reports [3] m d ^ olc = 1.35 ± 0.17 GeV, and 
other recent results from the K2K SciFi [TJ, K2K SciBar [7j and MINOS [S] collaborations 
similarly find central values higher than the above-mentioned world average. Quasielastic 
neutrino-nucleon scattering ((T]) is a basic signal process in neutrino oscillation studies. It 
is essential to obtain consistency between experiments utilizing different beam energies, and 
different nuclear targets. 

While a number of effects could be causing this tension, we here investigate perhaps the 
simplest possibility: that the parameterizations of the axial- vector form factor in common use 
are overly constrained. Such a possibility seems natural, considering that the dipole ansatz 
has been found to conflict with electron scattering data for the vector form factors. We do 
not offer new insight on whether other effects, such as nuclear modeling, could also be biasing 
measurements. However, we point out that by gaining firm control over the nucleon-level 
amplitude, such nuclear physics effects can be robustly isolated. 

The axial mass parameter as introduced in ([2|) is not well-defined, since the true form factor 
of the proton does not have a pure dipole behavior. Sufficiently precise measurements forced to 
fit this functional form will necessarily find different values for m^ ipole resulting from sensitivity 
to different ranges of q 2 . Let us define the axial mass parameter in terms of the form factor 
slope at q 2 = 0: = [F^(Q) /2Fa(0)}~ 1 ^ 2 . This definition is model-independent, and allows 
us to sensibly address tensions between different measurements. To avoid confusion, whenever 
([2|) is used we refer to the extracted parameter as m A ipole . We will show that the slope at 

1 The difference 0.055 is a correction to the conventionai representation of the pion electroproduction 
amplitude, as predicted by heavy baryon chiral perturbation theory [6]. 
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q 2 = is essentially the only relevant shape parameter for current data at Q 2 < 1 GeV 2 , and 
introduce the formalism to systematically account for the impact of other poorly constrained 
shape parameters on the determination of m^. A related study of the vector form factors of 
the nucleon was presented in [9]. 

The paper is structured as follows. In Section [2] we discuss the application of analyticity and 
dispersion relations to the axial- vector form factor of the nucleon. Section [3] presents results 
for the extraction of the axial- vector form factor slope from MiniBooNE data. We illustrate 
constraints imposed by our analysis on nuclear models, by determining the binding energy 
parameter in the Relativistic Fermi Gas (RFG) model of Smith and Moniz [IE] . Section [4] gives 
an illustrative analysis of constraints on the axial mass parameter from pion electroproduction 
data. Section |5]discusses the implications of our results. For completeness, Appendix lAl collects 
formulas for the RFG nuclear model. 

2 Analyticity constraints 

This section provides form factor definitions and details of the model-independent parameter- 
ization based on analyticity. 

2.1 Form factor definitions 

The nucleon matrix element of the Standard Model weak charged current is 



+ ri,F A (q 2 ) + -L q ^F P (q 2 ))u( n \p) , (3) 
m N J 

where g M = p' 11 — and we have enforced time-reversal invariance and neglected isospin- 
violating effects as discussed in Appendix |A] The vector form factors Fi(q 2 ) and F 2 (q 2 ) can be 
related via isospin symmetry to the electromagnetic form factors measured in electron-nucleon 
scattering. At low energy, the form factors are normalized as Fi(0) = 1, F 2 (0) = fi p — p n — 1. 
For definiteness we take a common nucleon mass, = (m p + m n )/2. Parameter values 
used in the numerical analysis are listed in Table [21 In applications to quasielastic electron- or 
muon-neutrino scattering, the impact of Fp is suppressed by powers of the small lepton-nucleon 
mass ratio. For our purposes, the pion pole approximation is sufficient Jl 

F P (q 2 ) « -^- 2 F A (q 2 ) . (4) 

The axial- vector form factor is normalized at q 2 = by neutron beta decay (see Table [2]). 
Our main focus is on determining the q 2 dependence of FA{q 2 ) in the physical region of 



2 Here and throughout, — 140 MeV denotes the pion mass. 
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quasielastic neutrino scattering, Q 2 = —q 2 > 0. As discussed in the Introduction, an expansion 
at q 2 = defines an "axial mass parameter" rriA, via 



-T<f + ■■■ 



m 



F A (q 2 ) = F A (0) 
Equivalently, we may define an "axial radius" r A , via 

F A {q 2 ) = F A (0) \l + r \> 

o 



m A 



2F A (Q) 
F A (0) ' 



Fa(0) 



(5) 



(6) 



The factors appearing in (JS]) and are purely conventional, motivated by the dipole ansatz 
(J2J), and by the analogous charge-radius definition for the vector form factors. Asymptotically, 
perturbative QCD predicts [TTJ1, [Tl] a ~ 1/Q 4 scaling, up to logarithms, for the axial- vector 
form factor. However, the region Q 2 < 1 GeV 2 is far from asymptotic, and the functional 
dependence of F A (q 2 ) remains poorly constrained at accessible neutrino energies. 



2.2 Analyticity 
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Figure 1: Conformal mapping of the cut plane to the unit circle. 



We proceed along lines similar to the vector form factor analysis in [9]. Recall the dispersion 
relation for the form factor, 



F A (t) 



to 



°° dt , lmF A (t' + iO) 
t'-t 



(7) 



where t = q 2 and the integral starts at the three-pion cut, t cut = 9m 2 ^. We can make use 
of this model-independent knowledge by noticing that the separation between the singular 
region, t > t cut , and the kinematically allowed physical region, t < 0, implies the existence of 
a small expansion parameter, \z\ < 1. As illustrated in Fig. HJ by a standard transformation, 
we map the domain of analyticity onto the unit circle in such a way that the physical region 
is mapped onto an interval: 



z[t, t cut , to) 
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t = 


t = t° pt (1.0GeV 2 ) 


\\F A \\2/\F A (t )\ 


1.5-1.7 


1.9-2.3 


\\F A \\oo/\F A (t )\ 


1.0-1.4 


1.4-1.8 



Table 1: Typical bounds on the coefficient ratios y/Ylk a k/ a o (first line of table) and \dk/ao\ 
(second line) in an axial- vector dominance ansatz. The range corresponds to the range 250 — 
600 MeV for the a\ width and the range 1190 — 1270 MeV for the a\ mass. 



where to is a f ree parameter representing the point mapping onto z = 0. Analyticity implies 
that the form factor can be expressed as a power series in the new variable, 

oo 

F A (q 2 ) = J2^z(q 2 ) k . (9) 

k=0 

The coefficients are bounded in size, guaranteeing convergence of the series. Knowledge of 
Im Fa over the cut translates into information about the coefficients in the z expansion [9] . In 
particular we have 

a = - [ d6Re F A [t(6) + zO] = F A (t Q ) , 

7T 







a fc >! = ~ [* d9lmF A [t(9)+i0} sm{k6) = \ [°° -JO-J *™* - to lmF A {t) sm[k9(t)} , (10) 
where 



TT Jo n Jt cut t-t \ t- t cut 



t = t + 2 f cut ^ = t(g). (11) 

1 — COS V 

2.3 Coefficient bounds 

For a given kinematic range < — t < Q^^, we can choose the free parameter to i n 
(JHJ) to minimize the resulting maximum size of \z\. It is straightforward to see that the 

"optimal" value of t is tQ Pt = t cut ^1 — + Qmax/^cut j , an d for this value of t , \z\ < 

K 1 + QLxAcut) 174 - !]/[(! + QLxAcut) 174 + 1] • For example, if the kinematic range is 
Qmax ~ lGeV 2 , then our expansion parameter is constrained to be \z\ < 0.2. Terms be- 
yond linear order in the expansion are suppressed by \z\ 2 < 0.04, etc., and are not tightly 
constrained by current experimental data. This is the sense in which the slope of the form 
factor (conventionally taken at q 2 = 0) is essentially the only relevant shape parameter. The 
effects of the higher order terms must of course be accounted for in assessing the uncertainty 
on extracted observables. We now turn to this question. 

The expansion coefficients appearing in (Q can be used to define norms, 

\\f a \\ p = • ( 12 ) 



In particular, \\Fa\ 



sup fc |a fc | 



lim 



| -Pa Hp provides a bound on the maximum coeffi- 



cient size. The finiteness of the integral appearing in the relation 
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together with ||-Fa||oo < H-^klb, establishes that a finite upper bound exists for the coef- 
ficients. As a first approach to estimating the actual bound ||Fa||oo> consider an "axial- 



vector dominance" ansatz, Fa ~ m ai /i. m a 1 — t — iT ai m ai ), where 



■in 



1230(40) MeV and 



250 — 600 MeV are the mass and width of the lowest lying axial- vector, iso- vector 



meson [12]. More precisely, let us define the form factor via its dispersion relation with [13 

A/v<r ai 



lmF A (t + iQ) 



(t - m 2 ai ) 2 + T 2 ai m 2 ai 



6{t - 1, 



cut J J 
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where M is a normalization constant determined below. Using the dispersion relation ([7]) with 
dUD we find, 



F A (t + i0) ■ 

where b(t) = t 
the values for 
show that 



M mlT ai 
n \b(t)\ 2 



llo g f l^)l 2 2 

2 V tr.iit. — t\ 



+ 



in 



ai 



t 



arg[fe(t cut )] + mB{t - t c 



(15) 



m 2 ai + iF 
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■ ai , and M is determined by the value of -Fa(O). Table CD displays 
and 1 1 -Fa | |oo computed in this ansatz. For the latter quantity one can 
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While this model is not a rigorous description of the true spectral function in (j7J), it indicates 
an order unity bound on the coefficients appearing in (Q. Additional support for an order 
unity bound is provided by a related detailed study of nucleon vector form factors [9], and by 
form factor studies in a wide range of meson transitions [T4"l [T5j . 

In the following numerical analysis, we follow [H] , and investigate fits with various bounds 
on coefficients, e.g. |afe| < 5 and \a k \ < 10. 



3 Extraction of the axial mass parameter 

The MiniBooNE collaboration has presented binned results representing the double differential 
cross section, da /dE^d cos 9^, for the quasielastic scattering process JI} on a neutron bound 
inside 12 C. We apply our description of F^(g 2 ) to extract (equivalently, ta) from the 
neutrino scattering data, under the assumption of a definite nuclear model, the Relativistic 
Fermi Gas model [16] as described in Appendix |A] 

Our theory prediction is obtained using f )54|) . integrating over the energy-dependent flux 
from Table V of [3]; this result is divided by 6 to obtain the per-neutron event rate, and divided 
by the total flux to obtain the flux-averaged cross section. Corresponding experimental values 
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Parameter 


Value 


Reference 


\v ud \ 


0.9742 


m 


Hp 


2.793 


m 




-1.913 


m 




0.1057 GeV 


m 


G F 


1.166 x 10~ 5 GeV" 2 


m 


m N 


0.9389 GeV 


m 


Fa(0) 


-1.269 


m 


£b 


0.025 GeV 




Pf 


0.220 GeV 





Table 2: Numerical values for input parameters. 



for the double differential cross section are taken from Table VI of [3J. We form an error 
matrix, 

Eij = (SatfSy + (5N) 2 <Ti(Tj , (17) 

where <7j = (da / dE '^d cos 9 ^AE^ A cos# M denotes a partial cross section, 5ai denotes the shape 
uncertainty from Table VII of |3j, and SN = 0.107 is the normalization error from [3j. We 
form the chi-squared function 

X 2 = 5>r Pt ' - eT cmy )E^(af^ - af cory ) , (18) 

ij 

and minimize \ 2 to find best fit values for m^. Error intervals are defined by A% 2 = 1. 
The nucleon form factors and the nuclear model employ parameter values listed in Table [2J 
Following the analysis of |3J, the vector form factors F\ and F 2 are given by the BBA2003 
parameterization [18]. We use a default value e& = 0.025 GeV, as extracted from electron 
scattering data on nuclei in |T7]. This value is different from the central value adopted in the 
MiniBooNE analysis [3], where e& = 0.034 ± 0.09 GeV. We show below that such a high value 
of €b is not favored by the MiniBooNE data, but investigate fit results for different values of 

The slope at q 2 = 0, and hence rriA from (jSJ) is most sensitive to low-Q 2 data. We analyze 
this sensitivity by considering the effect of a cut on Q 2 . The value of Q 2 for a given value of 
the observed muon energy and angle can be reconstructed assuming quasielastic scattering on 
a free neutron, but is not determined unambiguously once nuclear effects are included. As a 
proxy for Q 2 , we define an approximate "reconstructed" Q 2 , 



QL = 1K C E» - 2E?>y/E* - ml cos 6, - m 2 , (19) 
where E v u cc approximates the neutrino energy in the nucleon rest frame, 

m N - E/j,+ y/E^ - m 2 cos 9^ 
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We note that Q 2 CC coincides with Q 2 ec used by K2K in the limit e& — > [T] , and with Qq E used 
by MiniBooNE in the limit — )■ and equal proton and neutron masses [3J. For simplicity 
we have chosen to make the cut independent of the binding energy used in the nuclear model. 
We emphasize that this choice is used simply to define the subset of data to be analyzed, and 
does introduce theoretical uncertainty in the numerical results. 
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Figure 2: Extracted value of versus Q max . Dipole model results for m A ipolc are shown by 
the red circles; z expansion results with \a,k\ < 5 are shown by the blue squares, z expansion 
results with \af.\ < 10 are shown by the green diamonds. 



Our results are displayed in Fig. [21 where we compare extractions of mj^ e in the dipole 
ansatz (j2J) with extractions of uia employing the z expansion ([9]). We present results for data 
with Q 2 CC < g max , where Q 2 CC is defined in (02]) and Q max = 0.1, 0.2, . . . , 1.0 GeV 2 . We study 
two different coefficient bounds, |afc| < 5 and \a,k\ < 10. For definiteness we have truncated 
the sum in ([9]) at fc max = 7, but have checked that the results do not change significantly 
if higher orders are included. As Fig. [2] illustrates, the z expansion results lie systematically 
below results assuming the dipole ansatz. In contrast to results from the one-parameter dipole 
ansatz, high-Q 2 data have relatively small impact on the model-independent determination of 
rriA- Taking for definiteness Q^s* — 1-0 GeV 2 , we find 

m A = 0.85+^q 2 ± 0.09 GeV (neutrino scattering), (21) 

where the first error is experimental, using the fit with \a^\ < 5, and the second error represents 
residual form factor shape uncertainty, taken as the maximum change of the la interval when 
the bound is increased to |ajd < 10. As a comparison, a fit assuming the dipole form factor, 
and the same Q 2 max yields m d A ipolc = 1.29 ± 0.05 GeV| 

It is not our purpose in this paper to investigate in detail the additional uncertainty that 
should be assigned to (I2ip due to nuclear effects. We note that a fit of the MiniBooNE data to 

3 A dipole fit including the entire dataset without a cut on Q^ cc yields m A ipolc = 1.281q'q4. 
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the RFG model with free parameter e& yields the value, without an assumption on the value 
of m A , (for Q 2 max = 1.0 GeV 2 , k max = 7) 

e b = 28 ± 3 MeV , (22) 

where the result is insensitive to the choice of bound, \a k \ < 5 or \a k \ < 100 While the data 
do not appear to favor significantly higher values of e&, we note that for e& = 34 MeV [3], the 
result d2T]) becomes m A (e b = 34 MeV) = 1.05+°;^ ± 0.12, compared to m^ ipolc (e b = 34 MeV) = 
1.44 ±0.05. 

We have performed fits at different values of the parameter to, finding no significant devia- 
tion in the results. The results do not depend strongly on the precise value of the bound (e.g. 
|dfc| < 5 versus |ajt| < 10). Similar to [9], we conclude that the estimation of shape uncer- 
tainty in 

d2D» should be conservative. The fit (JUJ yields coefficient^ a = F A (0) = -1.269, 
d\ = 2.91}'q, a2 = —8^3. These values are in accordance with our assumption of order-unity 
coefficient bounds. As discussed in the Introduction, current experiments do not significantly 
constrain shape parameters beyond the linear term, a\. 

1.4 
1.2 
L.0 

-F A (-Q 2 ) o.S 

0.6 
0.4 

0.2 

0.0 0.2 0.4 0.6 0.8 1.0 

Q 2 

Figure 3: Comparison of the axial-vector form factor F A as extracted using the z expansion 
(green diamonds) and dipole ansatz (red circles). 

Figure [3] compares the form factor extraction resulting from the z expansion fit to the 
extraction from the dipole fit. Here we take Q max = 1.0 GeV 2 , /c max = 7 and \a k \ < 10 for the 
z fit. The dipole fit assumes m^ ipolc = 1.29 ± 0.05 GeV. 

4 Comparison to charged pion electroproduction 

The axial- vector component of the weak current defining F A (q 2 ) in can also be probed in 
pion electroproduction measurements. The electric dipole amplitude for threshold charged- 

4 Using a dipole ansatz for Q^ ax = 1.0 GeV 2 without fixing m^ lpol ° yields ^ = 22 ± 7 MeV. 

5 For this purpose we take A: m ax = 7 in ^ and enforce \ak\ < 10 for k > 3. 
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Kharkov (1978) [24] ' F5=j ■ ' 

Darosbury (1975/1976) [22] [23] 1 ■ i=i=i ' 

DESY (1973) Hi] 1 i * ' 

Frascati (1972) |20] 1 ■ I— a — i 

Average hn 

0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

m A (GeV) 

Figure 4: Extraction of tua using charged pion electroproduction measurements, in the dipole 
ansatz and in the z expansion. Datasets are as described in the text. Dipole results are shown 
as the red circles, and z expansion results with \a k \ < 5 are shown as the blue squares. 



pion electroproduction obeys a low-energy theorem in the chiral limit relating this amplitude 
to the axial- vector form factor of the nucleon [TS]. After applying chiral corrections, such 
measurements can thus in principle be used to determine tua- Data for this process have 
been interpreted in the context of the dipole ansatz (J2J). We found that the dipole assumption 
can strongly bias extractions of rriA in neutrino scattering measurements. In order to gauge 
whether the same statement is true for the electroproduction data, let us apply the z expansion 
to extract rriA from the inferred Fa{q 2 ) values for an illustrative dataset, taken from Refs. [201 
|2TJ [22j [231 121] • We have selected datasets that appear in the compilation [6] (cf. Figure 1 of that 
reference), and that also explicitly list inferred values of FA^q 2 ) (see also [251 I2S1 [271 123 129]). 
Figure [4] displays extractions of tua in both the z expansion and the dipole ansatz ([2]) for each 
of the five datasetsj§ For the larger bound \a k \ < 10, the slope of Fa(q 2 ) is not constrained to 
be positive by each individual dataset, and we display only the result for \a k \ < 5. Applying 
the z expansion to the entire (17 point) dataset, we find 

m A = 0.92+°;}^ ± 0.08 GeV (electroproduction) , (23) 

where the errors are experimental, and from residual shape uncertainty, as in f[2~T]) . In contrast, 
a fit of the same data to the dipole ansatz yields m^ ipolc = 1.00 ± 0.02 GeV. These averages 
are also displayed in the figure. We emphasize that our chosen dataset is not exhaustive We 
have not attempted to address questions such as correlations between different datasets, or 

6 For dcfinitcness, where necessary we have chosen one amongst different modefs for applied hard-pion cor- 
rections: the BNR prescription [30] in [22j [23j EI] , and the BNR prescription with first form factor assumption 
in [20] ("i^ = Ff" in Table 2 of EH] ). We have combined the low-Q 2 and high-Q 2 data from [22] and [23] to 
obtain the Daresbury(1975/1976) data point in Fig.[U 
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uncertainties from mo del- dependent hard-pion corrections. We leave a more detailed treatment 
to future work. 

5 Summary 

We have presented a model independent description of the axial-vector form factor of the 
nucleon. This form factor plays a crucial role in neutrino quasielastic scattering at accelerator 
energies, which is a basic signal process for neutrino oscillation studies, and is an important 
ingredient in normalizing the neutrino flux at detector locations. Recent tensions between 
measurements in neutrino scattering at different energies, and between neutrino scattering 
and pion electroproduction measurements indicate a problem in our understanding of this 
elementary process. 

Several studies have tried to address these discrepancies. Modified nuclear models [3TJ [32j 
[33] have been used to find an axial mass close to the MiniBooNE result. Other nuclear models 
include effects of multi- nucleon emission [MJ |35j [36j E7J [381 EH], and have been reported 
to obtain better agreement with the differential MiniBooNE data from [3J. One of these 
studies [3H] reports a dipole axial mass extracted from MiniBooNE data in agreement with 
world averages from [6J, [5]. Another group [ID], modifies the magnetic form factor Gm for 
nucleons bound in carbon but does not change the form factors Ge or Fa- The assumption of 
the dipole ansatz (j2J) is a crucial element in many of these studies^ Our analysis shows that 
this ansatz introduces a strong bias in measurements, which must be addressed in order to 
disentangle nucleon-level interactions from nuclear effects. 

Under the assumption of a definite nuclear model (the RFG model, summarized in Ap- 
pendix [Aj with parameter values as in Table [2]), we extract as defined model- independently 
in ([5]) from the differential MiniBooNE data [3J. The result is displayed in (1211) . = 
0-85lo 07 i 0-09 GeV. This result may be contrasted with a fit to an illustrative dataset for 
pion electroproduction displayed in (123"]) . m A = 0.921q;{3 ± 0.08 GeV. These values may be 
compared to fits using the dipole ansatz ([2]): m^ ipolc = 1.29 ± 0.05 GeV (neutrino scattering) 
and m^ ipolc = 1.00 ± 0.02 GeV (electroproduction). A discrepancy is apparent in the dipole 
ansatz (T5]), but can be ascribed to the unjustified and restrictive assumption on the form fac- 
tor shape. After gaining firm control over the nucleon-level amplitude, nuclear effects can be 
robustly isolated. For example, in the context of the RFG model, we extract the result ( 1221) 
for the binding energy parameter e&. 

The axial mass parameter, or equivalently, the axial radius (JDJ, is a fundamental parameter 
of nucleon structure. The results ([2~T]) . ([2"3]) can be expressed as 

J0.80lo;?7±0.12fm (neutrino scattering) 
TA ~ \ 0.74±g;J| ±0.05 fm (electroproduction) 

More precise measurements in both neutrino scattering and pion electroproduction are nec- 
essary to substantially reduce the errors on m^, or equivalently r^. This would be necessary 

7 A parameterization that modifies the dipole behavior at large Q 2 is presented in [41] . 
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to provide a model-independent confirmation of the convergence of chiral perturbation theory 
corrections based on comparison of electroproduction and neutrino scattering data. 

A related study of the nucleon vector form factors was presented in [9]. As described 
there, different expansion "schemes" are possible. For example, we may replace (j9j with 
(p{t)FA{t) = J2k akZ (t) k ' where (ft is analytic below t cut . A choice such as ~ (1 — t/m' 2 ) n 
with vn! ~ GeV could be used to enforce a 1/Q 2n falloff for asymptotic Q 2 , while retaining the 
known analytic structure of the form factor. Such modifications do not significantly impact 
the extraction of tua, arid we have focused on the simplest choice (to = and 0=1). 

Our study indicates that the error on the axial mass parameter extracted using the dipole 
ansatz is underestimated. While the errors from a model-independent analysis may be larger, 
it is essential to study model-independent numbers in order to draw firm conclusions. The 
simulation of more complicated neutrino scattering processes (e.g. pion and photon produc- 
tion), is indirectly affected by enforcing agreement with the quasielastic data. It is important 
for current and future neutrino experiments [121 H31 HU H21 El HH1 [31 H3 HH 09] to converge on 
consistent values for fundamental neutrino cross sections. 

The analysis presented here can be applied to other neutrino scattering datasets, involving 
different nuclear targets, and including neutral current scattering and antineutrino scattering. 
It is interesting to extend the analysis of electroproduction data; more precise low-energy 
electroproduction measurements have potential to impact the interpretation of future neutrino 
measurements. It is also of interest to incorporate model-independent constraints into more 
sophisticated nuclear models. 
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A Appendix: RFG model for quasielastic neutrino scat- 
tering 

A number of notations and conventions for the form factors and RFG nuclear model [TB] exist 
in the literature. For completeness we collect here the relevant formulas used in our analysis. 

A.l Nucleon matrix element of the weak current 

The relevant part of the weak-interaction Lagrangian is 




G F 



(25) 



The cross section for u(k) + n(p) — > £ (k') + p(p') on a free neutron is 



(26) 
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where the spin-averaged, squared amplitude is 

W\ = £ - l 5 )d\n(p))(p(p')\u lu (l - l5 )d\n(p))*. (27) 

spins 

The leptonic tensor neglecting the neutrino mass is (e 0123 = — 1) 

IT = 8(k"k' v + k v k'" - g^k ■ k' - ie^ pCT k p k' a ) • (28) 
The hadronic matrix element appearing in (127|) is parameterized by 

(p(pOI*fr(l - 7fe)d|r»(p)> = « (p) (p')r,(9)« W b) , (29) 
where q = k — k' = p' — p and we have defined the vertex function 

T,(q) = 7 ,Fi(g 2 ) + ~^a^F 2 (q 2 ) + -^F s (q 2 ) + ^ l5 F A (q 2 ) + ^±^ l5 F T (q 2 ) 
2m N m N m N 

+ ^ 5 F P (q 2 ). (30) 
m N 



We may write the cross section of (126]) as 

-^SS/w 2 ^' (31) 

where the nucleon structure function is 

<Wp^W« w + *)*'''■ (32) 

The hadronic tensor is 

= Tr[(f + m p )T^q){i) + m n )f ,(g)] , (33) 

where as usual, T = 7°r^7°. We may similarly analyze antineutrino scattering, u(k) +p{p) — > 
£ + (k ; ) + n(p'), using (l3Tj) . taking — > L^, and making the replacements m n -h- m p , r M (g) — )■ 
r M (-g) in H^. 

Imposing time- reversal invariance shows that Fi(q 2 ) are real. We will assume isospin sym- 
metry in the following, in which case F$ and Ft vanish, m n = m p = m^, and T^(—q) = T^q). 
The hadronic tensor has the time-reversal invariant decomposition 

ti^v — —9u,vti\ H — ^ — n~P q -"3 H 2~-"4 H ^ — n -"5 • 1-j4) 

2m; my 2m z N 

The ifj's are expressed in terms of the form factors Fi as 
H x = 8m 2 N F 2 - 2q 2 [(F, + F 2 ) 2 + F 2 } , 

H 2 = H 5 = 8m 2 N (F 2 + F 2 ) - 2q 2 F 2 , 

H 3 = -16m 2 N F A (F 1 + F 2 ) , 

H, = (F 2 + AF 2 ) - 2m 2 N F 2 - Am 2 N {F X F 2 + 2F A F P ) . (35) 
Expressions for complex Fi and nonzero F$,Ft can be found, for example, in [50] . 
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A. 2 Model for the nuclear matrix element 



We employ a standard treatment of nuclear effects, the "Relativistic Fermi Gas" (RFG) model 
as presented by Smith and Moniz in [16], based on the model presented in [51]. 

We assume that there are A nucleons inside the nucleus, with A/2 neutrons and A/2 
protons. The incoming neutrino interacts with a neutron with 3-momentum p, determined by 
some distribution rii(p). The final state proton phase space is limited by a factor of [1 — rif(p')] 
enforcing Fermi statistics. Symbolically, 

Onuclcar = <8> Ofree(P P') ® [1 ~ n f {p% (36) 

and more explicitly 

/d^p 

{i^/(^/(^( 2 ^-^' + «^}t 1 -"'Wi' < 37 > 

To arrive at the final model, two modifications are made. First, we make the replacement 
k -p — > Ej-Ep in the prefactor of ( |37j) . This replacement ignores a correction from the nonzero 
velocity of the initial state nucleon. It corresponds to the model of [IB] , adopted by [3]; for 
definiteness we have followed this convention. Second, we incorporate a "binding energy", e&, 
by expressing H^ u as a function of Lorentz 4- vectors p M , g M as in (1341) and then making in fl37|) 
the replacements 

p° -> e p = E p - e b , p'° -> e p , = £ p , , (38) 

with £"p = \Jm 2 N + |p| 2 . Again, there is some arbitrariness to the insertion of e b into the 
formalism; for definiteness we have followed the conventions of HBJ. The cross section is then 



^nuclear — j / , xo oz7 , , (39) 

l6\k-p T \ J {2ny i 2E kl 

where pj, is the 4-momentum of the target nucleus with mass mx = ArriN(l — e&). We work in 
the target rest frame where pj. = ttitSq. The model nuclear structure function is defined 

as 

W, u = Jd 3 pf(p,q ,q)H, u (e p ,p-q°,q), (40) 

with 

f{P,q ,q) = -7-2-rn(p)[l-n f (p + q)] . (41) 

^ t7r tptp_|_q, 

The distribution of neutrons and protons is 

rn(p) = 0(Pf ~ \p\) , n f (p') = 6(p F -\p'\), (42) 
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where pp is a parameter of the model. The normalization V is fixed by requiring A/ 2 neutrons 
below the Fermi surface (accounting for 2 fermionic spin states), 



.4 



2V 



d 3 p 



l n AP) 



2 J (2tt^ 3 

We can expand W^ v in a similar way to H^ v in ([3 



V 



3n 2 A 
~2pJ 



(43) 



W, 



T T 
PJtPt 



-g„Wi + - i e ^f p ^w 3 + ^fw, + ^-f^W 5 . (44) 



2m 2 n 



The functions Wi are related to integrals over Hi. The relations can be expressed in terms of 
the following integrals [16]: 



ai= d 3 pf(p,q) 



P\ 



o 3 



d 3 pf(p,q) 



(P 



z\2 



m 



2 ' 
N 



a 2 = j d 3 pf(p,q) 

a 4 = / d 3 pf(p,q) -f 
J m ]\ 



2 ' 
N 



a 5 = / d 3 pf(p,q) 



£ P P 



m 



2 ' 
JV 



ae = / d 3 pf(p,q) 



p 

m N 



a 7 



d 3 pf{p,q) — , 

m N 



(45) 



where \p\ 2 = (p x ) 2 + (p y ) 2 + (p 2 ) 2 and the z axis is parallel to q. A straightforward but tedious 
comparison shows that 



Wi = a x Hx + -(a 2 - a 3 )H 2 



W 2 

w 5 



u lu 1 , 

\q\ 2 \q\ |ql 2 

— a 7 _ i — T a 6 #3 , 

mjv v \q\ 



[a 2 - a 3 ) 



#2 



rrin 



m 



N 



a-i_H A + ——a 6 H 5 + — ^ (3a 3 - a 2 )# 2 



°7 _ i — r a 6 LH5 + -j— r 

v \q\ J \q\ 



2a 5 + ^-Aa 2 - 3a 3 ) 
9 



#9 



(46) 



where we are using u = q°. Recall that the Hi are functions of q 2 = u 2 — \ q\ 2 . For the integrals 
ai let us define = oj — e^, and observe that 



S(e p - e p+q + q°) = 5(E P - E p+q + u, 



E. 



cff 



\P\\Q\ 



^5 cos9 



'pq 





<? 


2 


4- 2u e gE p 


2 


P\ 


\q 





• (47) 
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The integrals dj can be expressed in terms of 



27i\q\ 



E p — €b V m N 



E r 



(4£ 



for j — 0, 1, 2. In particular 



(£ + e 6 log(£-e 6 )) 
1 



^lo 



27rmAr|qr| 



E 2 + e b (E + e b \og(E - e b )) 



2iTm N \q\ [3 



-E 2 + e b (E + e b \og(E-e b )) 



(49) 



Up to an overall constant these are the fej's of [16]. Introducing c = —oj e g/\q\, d = — (u 2 s — 
\q\ 2 ) / {2\q\rriN) , we can express the Oj's as 



a x = b , a 2 = b 2 — b Q , a 3 = c 2 b 2 + 2cdb\ + d 2 & , 04 = ^2 



2e6 



&1 



3. 



fro 



a 5 = -cb 2 



-^—c — d ) bi + -^-db , ct6 = — — ^0 5 a-7 = b\ — &o 

772 TV / ty 



The range of integration is restricted by the conditions, 



E p < E F = Jm 2 N +p 2 F < E p+q = Ep + uj, 



eff 



1 < 



off 



|<j| 2 + 2u eS E f 



< 1 



2| g |«/£5-m; 



(50) 



(51) 



iV 



The latter condition can be expressed as 



E p cd+Vl 



d 2 



m N 



1-c 2 



E p 

m N 



cd - VI - c 2 + d 2 
1-c 2 



> 



Define 



£"10 = max ( E F — tu c g, tun 



cd+Vl-c 2 + d? 



Ehi = Ei 



(52) 



(53) 



1-c 2 

Then if E\ Q > E^, there is no contribution for the given kinematics. 

In the rest frame of the nucleus, let Eg and \Pe\ = \J E 2 — m 2 be the energy and 3- 
momentum of the charged lepton, and let 6e be the angle between the 3-momenta of the 
leptons. From ( )39l) . the final expression for the differential cross section of neutrino- nucleus 
scattering is 



da 



nuclear 



G%\Pf 



dEgd cos 0£ 1 Qir 2 rriT 



2{E e - \P t \ cos9 e ) W l + (E t + \P t \ cos 9 e )W 2 
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± 



in 



ml, 



(E e - \P e \ co&O t ){E v + E t ) - mj W 3 + -L(E e - \P t \ cos6 e )W 4 - -L W 5 }, (54) 



in-. 



where W\ are given in (1461) . and where the upper (lower) sign is for neutrino (anti- neutrino) 
scattering. 
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